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ABSTRACT 



Three models of an aircraft repair facility are developed 
for use in computing manpower savings achievable by scaling- 
up aircraft maintenance shops. One model is a bivariate 
Markov process model requiring a Gauss-Seidel iterative 
algorithm for its solution. The other two are the simple 
Repairman problem and the M/M/S queue, the solutions to which 
can be found in most introductory texts on stochastic models. 
The simple models are compared to the more refined bivariate . - 
model in regard to predictive accuracy. The simple models 
are used to compute the manpower savings achievable in one 
scaling-up scheme, and predict a 20% savings in maintenance 
pers onnel. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION 12 

A. BACKGROUND 12 

B. MOTIVATION 14 

C. USES 15 

1. Analysis of Current Structure 15 

2. Analysis of Alternative Structures — 15 

3. Application to Other Support 

Structures 15 

D. SCOPE 15 

II. DESCRIPTION OF CURRENT REPAIR ORGANIZATION — 18 

III. A BIVARIATE MODEL 25 

IV. TWO DECOMPOSITION MODELS 45 

V. COMPARISON OF MODELS 58 

VI. CALCULATION OF MANPOWER SAVINGS 92 

VII. CONCLUSIONS 96 

LIST OF REFERENCES 97 

INITIAL DISTRIBUTION LIST 98 



5 



LIST OF TABLES 



TABLE 

I. Number of Bivariate States vs. Number 

of Aircraft 28 

II. Data Conversion Table 59 

III. Measures of Congestion for Shop 1 Calculated 

from Three Different Models 62 

IV. Measures of Congestion for Shop 2 Calculated 

from Three Different Models 65 

V. Measures of Congestion for Shop 3 Calculated 

from Three Different Models 68 

VI. Measures of Congestion for Shop 4 Calculated 

from Three Different Models 71 

VII. Measures of Congestion for Shop 5 Calculated 

from Three Different Models 74 

VIII. Measures of Congestion for Shop 6 Calculated 

from Three Different Models 77 

IX. Measures of Congestion for Shop 7 Calculated 

from Three Different Models 80 

X. Measures of Congestion for Shop 8 Calculated 

from Three Different Models 83 

XI. Measures of Congestion for Shop 9 Calculated 

from Three Different Models 86 

XII. Measures of Congestion for Shop 10 Calculated 

from Three Different Models 89 

XIII. Manpower Savings and Measures of Congestion 

for a Consolidated Back-shop Facility 94 



6 



LIST OF FIGURES 



FIGURE 

1. Aircraft Maintenance Flow Diagram 19 

2. Set Element Interaction in Current Facility 24 

3. Set Element Interaction in Bivariate 

Model of Current Facility 26 

4. Transition Diagram of Hypothetical Shop 30 

5. Balance Equations for the Bivariate Model 31 

6. Set Element Interaction in Decomposition 

Model of Current Facility 4 7 



7. Transition Diagram of Hypothetical Repair 

Shop Modeled as Two Independent Finite 
Arrival Source Queues (Decomposition 



Repairman Model) 48 

8. Plot of Probability vs. Aircraft Down for 

Flight-line Repair at Shop Nr. 1 6 3 

9. Plot of Probability vs. Aircraft Down for 

Back-shop Repair at Shop Nr. 1 64 

10. Plot of Probability vs. Aircraft Down for 

Flight-line Repair at Shop Nr. 2 66 

11. Plot of Probability vs. Aircraft Down for 

Back-shop Repair at Shop Nr. 2 67 

12. Plot of Probability vs. Aircraft Down for 

Flight- line Repair at Shop Nr. 3 69 

13. Plot of Probability vs. Aircraft Down for 

Back-shop Repair at Shop Nr. 3 70 

14. Plot of Probability vs. Aircraft Down for 

Flight-line Repair at Shop Nr. 4 72 

15. Plot of Probability vs. Aircraft Down for 

Back-shop Repair at Shop Nr. 4 73 

16. Plot of Probability vs. Aircraft Down for 

Flight-line Repair at Shop Nr. 5 75 



7 



17. Plot of Probability vs. Aircraft Down for 

Back-shop Repair at Shop Nr. 5 

18. Plot of Probability vs. Aircraft Down for 

Flight-line Repair at Shop Nr. 6 

19. Plot of Probability vs. Aircraft Down for 

Back-shop Repair at Shop Nr. 6 

20. Plot of Probability vs. Aircraft Down for 

Flight-line Repair at Shop Nr. 7 

21. Plot of Probability vs. Aircraft Down for 

Back-shop Repair at Shop Nr. 7 

22. Plot of Probability vs. Aircraft Down for 

Flight-line Repair at Shop Nr. 8 

23. Plot of Probability vs. Aircraft Down for 

Back-shop Repair at Shop Nr. 8 

24. Plot of Probability vs. Aircraft Down for 

Flight-line Repair at Shop Nr. 9 

25. Plot of Probability vs. Aircraft Down for 

Back-shop Repair at Shop Nr. 9 

26. Plot of Probability vs. Aircraft Down for 

Flight-line Repair at Shop Nr. 10 

27. Plot of Probability vs. Aircraft Down for 

Back-shop Repair at Shop Nr. 10 



76 

78 

79 

81 

82 

84 

85 

87 

88 

90 

91 



TABLE OF SYMBOLS 



K 


total number of shops at the repair facility. 


M 


total number of aircraft at the repair facility. 


B 


as a subscript denotes back-shop activity. Also 
used as a random variable indicating the number 
of aircraft down at a shop for back-shop repair. 


F 


as a subscript denotes flight-line activity. 

Also used as a random variable indicating the 
number of aircraft down at a shop for flight-line 
repair. 


X 


arrival or failure rate 


R 


the set which characterizes the repair shops 


S 

s 


the number of repairmen at the shop or shop 
portion denoted by the subscript s 


V 


repair or service rate 


P 


the set of performance measures which characterize 
the repair facility 


E(arg) 


the expected value of the argument 


Var (arg) 


the variance of the argument 


N. 

1 


number of aircraft down for repair at the shop 
indicated by the subscript i 


N 

s 


number of possible states of the bivariate model 


Q i 


number of aircraft in the queue at the shop 
indicated by the subscript i 


w. 

i 


waiting time at the shop indicated by the 
subscript i 


D i 


delay time at the shop indicated by the subscript i 


p fb 


the long-run probability of exactly f aircraft 
down for flight-line repair and b aircraft down 
for back-shop repair 


f 


the number that the random variable F assumes 



9 



b the number that the random variable B assumes 

(k) 

P' the kth iteration value of P_. 

fb fh 

P^ the long-run probability of f aircraft down for 

flight-line repair 

P^ the long-run probability of b aircraft down for 

back-shop repair 

P the long-run probability of n aircraft down for 

repair (n = f + b) 

L the expected number of aircraft in the system 

the average arrival rate at the flight-line 

X_, the average arrival rate at the back-shop 

D3. 

X the average arrival rate into the shop, 

flight-line and back-shop both 



10 



ACKNOWLEDGMENT 



The author wishes to thank Professor Donald P. Gaver 
for his guidance and encouragement throughout the research 
and writing phases of this thesis. Thanks also go to 
Professor Richard W. Butterworth for his assistance as second 
reader. The following persons at the RAND Corporation 
provided data and information about aircraft logistics 
structures which is the core of the thesis: I. Cohen, 

S. Drezner, and C. Roach. The author is especially grate- 
ful to George W. Fromknecht, Captain, U.S. Navy (Retired) 
who is largely responsible for my interest in Operations 
Research and was instrumental in my selection for the 
Operations Research curriculum at the Naval Postgraduate 
School. 



11 



I. INTRODUCTION 



It is quite a three-pipe problem. 

-Sir Arthur Conan Doyle, The Red-headed League 

A. BACKGROUND 

Under the auspices of the United States Air Force Project 
Rand, the RAND Corporation is currently studying ways to 
improve the aircraft maintenance posture of the United States 
Air Force. In the past, RAND's Logistics Studies Program 
has been primarily concerned with the improvement of particu- 
lar logistics functions. They have since departed from this 
basically micro approach to logistics problems for several 
reasons. First, the Air Force has over the years developed 
an in-house capability to analyze micro-oriented logistics 
problems, so that RAND assistance in this area is no longer 
required to the extent that it was in the past. Second, 
with the Air Force logistics analysts now concentrating on 
the micro-oriented problems, RAND's analysts are free to 
concentrate their efforts on macro-oriented approaches to 
the total logistics system. "If the Air Force can take care 
of the trees, then let us study the forest," is a good, 
although casual, expression of the present RAND approach 
to logistics problems. This ivory tower attitude may result 
in the introduction of new logistics structures in the future 
that will dramatically reduce support costs. 
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This shift in emphasis from a micro approach to logis- 
tics function, to a concentration on macro-oriented problems, 
is a fairly recent development. In July, 1973 a project 
called ANALOGS 1980s (Analysis of Logistics Structures for 
the 1980s) was formally initiated as a part of RAND's 
Logistics Studies Program. Sponsored by Lieutenant General 
William W. Snavely, Deputy Chief of Staff, Systems and 
Logistics, Headquarters, USAF, the project has as its explicit 
goal the exploration of logistics structures that could 
dramatically reduce support costs without degrading mission 
capability. Later, on September 9, 1974, the Deputy Chief 
of Staff, Systems and Logistics, provided further impetus 
to the macro approach when he issued a Program Management 
Directive authorizing a Maintenance Posture Improvement 
Program (MPIP) with two broad objectives: 

1. To reduce costs of maintenance in peacetime and 
wartime contingencies . 

2. To increase the effectiveness of maintenance in 
support of operational missions. 

Information on ANALOGS 80s and on work to support the MPIP 
has been provided by S. Drezner and I. Cohen of RAND, to 
whom the author is grateful. 

As a result of RAND work done in support of MPIP and 
the ANALOGS 80s project, there is a current interest in 
two organizational structures which attempt to reduce support 
costs by scaling-up maintenance activities, both of which 
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represent a departure from the organizational structures 
presently employed by the USAF. The older structure of the 
two, called Queen Bee, specifies the centralization of cer- 
tain intermediate level aircraft maintenance at large central 
bases called Queen Bees. The second structure, called the 
Reallocation of Activities Alternative (RAA) , grew out of 
the ANALOGS 80s project analysis. RAA specifies the decen- 
tralization of combat-mission-oriented activities at Combat 
Mission Bases (CMBs) , and the centralization of logistics 
and training activities at Support Mission Bases (SMBs) . 

Costs and payoffs calculations for each of the two struc- 
tures are needed before a decision to implement is made 
regarding either structure. In late September, 1974, work 
began in earnest on the development of models which could 
be used to assess the effectiveness of such configurations. 

B. MOTIVATION 

At the outset of the modeling effort, it was felt that 
the models developed could be used in a broader context than 
the RAND study. The Air Force is not the only branch of the 
armed forces with an aircraft logistics problem; the Navy 
has one too. Additionally, there were similarities between 
the structures employed by the two branches, so that the 
models developed for the Air Force could be applied to Navy 
problems with little or no modification. There was, thus, 
strong motivation for a thesis focusing on the modeling of 
aircraft logistics structures. 
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C. USES 



1. Analysis of Current Structure 

Once models for the current structure have been 
developed, they could provide answers to several cogent 
questions. For example: Is the present manpower allocation 

within the current structure reasonably satisfactory, if 
not optimal. How might the manpower be reallocated to 
improve readiness? What degradation in readiness can be 
expected for various reductions in manpower? If manpower 
reductions are necessary, where will the effects of the 
cuts degrade readiness the least? 

Most importantly, models for the current structure 
also provide a basis for the comparison of the current 
structure with a proposed alternative structure. They can 
also be used to analyze the effects of increasing the number 
of aircraft that the current maintenance structure must 
support. Alternatively, they can be used to indicate 
acceptable manpower cuts in the event that the number of 
aircraft and/or the amount of flying is decreased. Provided 
that reliability estimates of a new proposed weapon system 
were available, the models would be useful in predicting 
manpower costs in support of a new weapons system. It may 
be that an alternative support system may be better suited 
for the proposed weapon system than the current support 
system. The models would help resolve problems such as these. 

2 . Analysis of Alternative Structures 

As stated previously, before decisions are made 
regarding alternative structures, the current system should 



15 



be modeled to provide a basis for comparison of the current 
structure with the proposed structure. In the RAND study, 
the current structure was modeled first. This provided 
measures of performance for the system. Then the alternative 
structure was modeled. The model of the alternative struc- 
ture was used to show that measures of performance on par 
with those attainable with the current structure were achieva- 
ble with the proposed structure with fewer personnel. 

There is a trap to watch for in this kind of com- 
parative analysis. If the manning of the current support 
system is not optimal, then the comparison with the alterna- 
tive structure, assumed operating optimally, injects a bias 
in favor of the alternative structure into the analysis. 

It may be that the reallocation of manpower within the current 
structure will yield measures of performance that the 
alternative structure cannot satisfy with less personnel. 

3. Application to Other Support Structures 

It was stated earlier that there were similarities 
between Air Force and Navy aircraft support structures. 

There are also similarities between aircraft support struc- 
tures and ship support structures. If there are economies 
attainable by scaling-up certain logistics functions within 
the aviation community, then similar economies may exist in 
the surface and submarine communities. 

D. SCOPE 

A logistics support structure has not been defined up 
to this point. A logistics support structure is defined as 
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the organization, training, quantification, and geographical 
positioning, where each is applicable, of military and 
civilian personnel, tools, equipment, armament, spare parts, 
fuel and physical plant facilities with the explicit pur- 
pose of maintaining a weapon system in a satisfactory state 
of readiness to perform peacetime missions and meet wartime 
or combat contingencies. This nebulous definition is pro- 
vided to point out the magnitude of the task confronting 
the analyst who aspires to model the structure iri to to . 

It was decided to look at a small but significant part 
of the total structure, and see if there were economies to 
be derived from the scaling-up scheme proposed in the Queen 
Bee and RAA structures . To restrict the analysis to the 
organization of repair personnel and calculate manpower 
savings achievable with the alternate structures appeared 
to be the best approach. This approach was appealing for 
several reasons. First, support personnel costs are a 
significant contributor to the total lifetime cost of a 
weapon system; a significant reduction here can dramatically 
reduce the cost of a weapon system. Second, repair organi- 
zations can be modelled as queueing systems. There are 
standard calculations of measures of congestion associated 
with queues which makes side-by-side comparison of two 
different repair organizations relatively simple. Third, 
many other costs associated with the logistics structure de- 
pend in large part on the numbers of repair personnel and on 
their organization and positioning within the logistics structure. 
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II. DESCRIPTION OF CURRENT REPAIR ORGANIZATION 



Why is this thus? What is the reason of this thusness? 

-Artemus Ward, Artemus Ward's Lecture 

The current aircraft repair organization, hereafter 
referred to as an aircraft repair facility, consists of K 
shops, designated Shop 1, Shop 2, and so on, up to Shop K. 
There are two basic categories of maintenance tasks performed 
at each shop: flight-line repairs and back-shop repairs. 

The flight-line repair category consists of all those repairs 
that can be made on board the aircraft at the flight-line. 

The back-shop repair category consists of all other types 
of repairs. The repair facility supports an operational 
unit which operates a fixed number of aircraft, equal to M. 

Figure 1 is a schematic diagram of the operational unit 
and a repair facility composed of six shops which shows how 
the facility and the operational unit interface. When a 
mission is scheduled, the required number of aircraft are 
removed from the pool and preflight maintenance is conducted. 
When preflight maintenance is satisfactorily completed, the 
aircraft conduct the mission, return to base where post- 
flight maintenance is conducted. Once this is completed, 
they are returned to the pool. During stand-down periods, 
other maintenance is conducted in addition to the preflight 
and postflight maintenance. This maintenance is referred 
to as scheduled maintenance. The portion of the diagrammed 
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repair shop with enclosed capital B represents the back-shop 
activities of the repair shop, and the portion with circled 
numeral and capital F indicates the flight-line activities. 
Small f and small b denote flight-line and back-shop failures. 
When a failure occurs, the cognizant shop dispatches a 
repairman to the aircraft where the failure occurred. If 
it can be repaired onboard the aircraft it is designated a 
flight-line failure. If not, then it is considered a back- 
shop failure. 

Back-shop or flight-line casualties can occur during 
missions, preflight, postflight or scheduled preventive 
maintenance. We assume that casualties do not occur to 
aircraft that are in the pool or already down for repair. 

It is also assumed that the tempo of flight operations is 
such that the number of aircraft in the pool is very small 
compared to the number of aircraft available to the opera- 
tional command, so that the total number of aircraft which 
can fail at any given time is approximately equal to the 
total number available, M, minus the number that are down 
for repair. 

A repair facility consisting of K shops, each capable 
of performing two categories of maintenance tasks implies 
that there are 2K possible failures that the facility can 
handle. It is assumed that the 2K possible failures exhaust 
all of the types of failures that can occur on the aircraft. 

It is further assumed that all aircraft at the operational 
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command exhibit identical failure characteristics; all 
aircraft fail into the flight-line portion of Shop 1 at one 
constant rate, and into the flight-line of Shop 2 at a 
different constant rate, and so on, for all K shops. Back- 
shop failures exhibit the same property. 

Each shop at the facility is characterized by the number 
of repairmen in the shop, and two repair rates, one for the 
flight-line repairs and the other for back-shop repairs. The 
repairmen are assumed to be equally capable at any task 
confronting the shop. This makes it possible to characterize 
the shop by specifying the number of repairmen, the flight- 
line repair rate, and the back-shop repair rate. Repairmen 
at one shop cannot assist repairmen at a different shop, so 
that if failures occur when all repairmen are busy at one 
shop a queue begins to form. 

A priority system within each shop is in effect to insure 
that repairmen are assigned to flight-line repairs before 
commencing work on back-shop repairs. Furthermore, every 
repair task can be accomplished by one and only one repairman. 

Let F denote the set of 2K different failure rates which 
characterize each of the M aircraft; then 



~ {X f1' X B1' X F2' X B2' 



*' X FK' X BK }f 



( 2 . 1 ) 



where the alphanumeric subscript indicates the type of 
failure, F and B denoting flight- line and back-shop 
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respectively, and the numeral denoting the cognizant repair 
shop number. 

Let R denote the set which characterizes the repair shops. 
Then, 

> ( 2 . 2 ) 

/ 

where is the number of repairmen, y F ^ is the flight-line 
repair rate, and y^ . is the back-shop repair rate at Shop i. 

Another set of interest contains measures of congestion 
at the shops in the facility; it is designated P and is 
defined as follows : 
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where 



P = {P i' P Fi ' P Bi } i = i' 2 ' 3 '---’* 



( 2 . 3 ) 
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The subscript i refers to the ith shop, and the subscripts 
Fi and Bi refer to the flight-line and back-shop activities 
of the ith shop. The variable N is the number in the system, 
and the variable Q is the length of the queue. W is the 
waiting time, and D is the delay time. Waiting time is 
defined as the total time the aircraft is down for repair. 
Delay time is the time spent waiting for repairs to begin. 

Figure 2 illustrates the interaction of the elements of 
set F with the set R, resulting in the performance measures 
in set P. Note the feedback loops out of R; these represent 
the interaction of the shops with each other. This inter- 
action is not well understood at present, and models of the 
present facility developed on the following pages ignore 
it. 
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SET ELEMENT INTERACTION 
IN CURRENT FACILITY 

FIGURE 2 





III. A BIVARIATE MODEL 



How happy could I be with either. 

Were t'other dear charmer awayl 

But while ye thus tease me together. 

To neither a word will I say. 

-John Gay, The Beggar's Opera 

A bivariate model is developed in order to obtain values 
for the measures of congestion contained in the set P. The 
basic approach in modeling the system with a bivariate 
model is to treat each shop independently, and model the 
facility as K different queueing systems, each queueing 
system maintaining M aircraft. Figure 3 illustrates the 
approach. 

Consider the queueing system defined by the interaction 
of Shop N with M aircraft. Assume that repair times and 
interarrival times are exponentially distributed. The state 
of the system at any time t can be represented by the ordered 
pair (F ( t) ,B(t)) , where F(t) is the number of aircraft down 
for flight-line repairs, and B(t) is the number of aircraft 
down for back-shop repairs at Shop N. If flight operations 
and maintenance are conducted throughout the 24 hour day, 
so that the process is not dependent on time, the bivariate 
model belongs to a class of probability models known as 
continuous -time Markov chains. These probability models 
are characterized by the fact that 

1. They can be completely described at any time by 
specifying their state at that time, and 
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SET ELEMENT INTERACTION IN 
BIVARIATE MODEL OF CURRENT FACILITY 

FIGURE 3 



2. the time until a transition is made from one state 



to a neighboring state is an exponentially distributed 



random variable. 



Let 



P fb (t) = P[F(t) = f, B ( t) = b] 



(3.1) 



Concentrating on long-run probabilities, of more interest 
is the quantity 



is the long-run or steady-state probability of exactly 
f aircraft down for repair at the flight-line, and b aircraft 
down for repair at the back-shop of the shop under considera- 
tion. Pf b is also the long-run proportion of time that there 
are exactly f and b aircraft down for repairs at the flight- 
line and back-shop respectively. 

For a queueing system consisting of one shop servicing 
M aircraft, the possible states of the system are (f ,b) , 
f = 1,2,3,..., M, b = 1,2,3,..., M, such that f + b <_ M. 

The probability distribution of interest is the steady-state 
probabilities associated with the possible states of the 
system 




(3.2) 



= lim P [F ( t) = f, B ( t) = b]. 



(3.3) 
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{P fb >, f = 1/2,3,... ,M, b = 1,2,3,... ,M, 
s ,t . f + b _< M . 

The number of possible states is a quadratic function of the 
number of aircraft serviced by the shop. Denoting the 
number of possible states by N , the relationship is 

. (M + 1) (M + 2) . . (3 . 4) 

s / 

The state space for a system with a large number of aircraft 
can be quite large, as Table I illustrates. 



TABLE I 

Number of Bivariate States 
Vs. Number of Aircraft 



M 


N s 


1 


3 


5 


21 


10 


66 


25 


351 


50 


1326 


75 


2926 


100 


5151 


200 


20301 



The general principle which enables one to determine 
the state probabilities is derived from a rate-equality 
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principle which applies to continuous -time Markov chains. 

The principle applied to the bivariate process states that 
for each possible state (f ,b) , the rate at which the process 
enters the state (f,b) equals the rate at which it leaves 
the state (f,b) . The application of this principle to each 
of the N possible states of the system will yield a set 
of balance equations (N in all) which must be solved for 

b 

the set of state probabilities 

Figure 4 is a transition diagram for a hypothetical 
shop consisting of three repairmen with repair responsi- 
bilities for five aircraft. This shop represents the repair 
facility described in Chapter II. The circles represent 
the possible states of the system. Arrows indicate the 
possible transitions between states, and the labels on the 
arrows represent the rates at which the process moves from 
one state to another. The transition diagram simplifies the 
task of setting up the balance equations. 

Balance Equations for a system of M aircraft failing 
into a shop with S repairmen are presented in Figure 5. 

These can be verified by referring to the transition diagram 
of Figure 4 . 

The solution of the system of balance equations can 
proceed algebraically utilizing a successive substitution 
scheme for small systems, but for large systems the solution 
becomes intractable. For this reason, an iteration algo- 
rithm was invoked so that f° r s Y steins °f practical 
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TRANSITION DIAGRAM OF HYPOTHETICAL SHOP 
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BALANCE EQUATIONS FOR THE BIVARIATE MODEL 
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FIGURE 



size could be found. The algorithm used was the Gauss- 
Seidel iterative method described by Varga in [Ref. 1] and 
by Young and Gregory in [Ref. 2], The application of itera- 
tive algorithms to Markov-type problems is an area of active 
research at the Naval Postgraduate School, as it is particu- 
larly well suited for the solution of large systems of 
balance equations . 

The algorithm as it applies to the problem at hand is 
set forth here without mathematical underpinnings. First, 
all balance equations are rearranged so that is expressed 

as a function of p (f _ 1)b / P (f+i)b' P f(b-1)' and P f (b+1) * 

This is easily done, as reference to equations (3.5) and 
(3.6) will show. Second, assume values for all P^* One 
scheme would be to assign P^ = 1/N g , for all f and b in 
the state space. Here, N g is the number of possible states. 
The closer the initial assignment of values of P^ is to the 
solution, the faster the algorithm will converge. Third, 
denote the assumed values of P^/ Pf b ^ • f° r t * ie zeroth 
iteration. Fourth, establish a scheme for incrementing sub- 
scripts f and b. The scheme used in the problem under con- 
sideration was to start at P QQ and increment f values by one 
until P^q was reached, then b was incremented by one, and 
f values were then incremented by one, starting at f = 0 
again. Thus, after P^q was reached, the incrementing scheme 
called for P^, p ^/ P 21' ••••' P (m- 1)1 ‘*' n t * lat or< 3er. The 
subscript b is again incremented by one and f values 
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incremented by one, so that f = 0,1,2,..., M-2 . This process 
is continued until all in the state space have been 

reached. This constitutes one iteration. So one iteration 
will consist of N g increments. Unless we do something at 
each increment, of course, the algorithm would be rather 
pointless, so on to the remaining steps of the algorithm. 

Denote the value of computed at the nth iteration 

P^k . Recall that all balance equations were rearranged 
so that P^ equals some function of p (£_!)]-,' P (f+l)b' 

P f(b-1) ' P f(b+1)’ Eac h balance equation will generate a 
different function, but to keep the notation simple, denote 
that function generated by the balance equation under con- 
sideration g, where it is understood that the gs will be 
different for each equation. We have 



as the nth iteration approximation of P.^* Note that it 
is computed from the (n-l)th iteration approximations of 



P fb g(P (f-l)b' P (f+l)b' P f (b-1) ' P f (b+1) 



(3.7) 



Define 




(3.8) 



P 



(f-l)b' 



P 



(f+l)b 



P 



f (b-1) 



and P 



f (b+1) * 
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The nth iteration values of are computed as follows: 

1. For the first in the incrementing scheme, compute 

P^ n ^ in accordance with equation (3.8) . 

2. Immediately change P^ 11 ^ to P^ 1 ^ • For example, if 

Pqq"^ is being computed, then immediately after the step 
where a value for p gQ^ is obtained, change p qq^ to that 
value. Now, when the incrementing scheme calls for the 
computation of a P^^ involving a P / the algorithm will 
use the most recent approximation of P Q0 , which is the 
value computed for • 

3. Continue through the incrementing scheme, doing steps 
1 and 2, until the algorithm has incremented through all 
possible states. 

4. Make the following computation: 



S 



M' M-b 
Z Z 
b=0 f=0 



(n) 

fb 



(3.9) 



( n ) 

5. Normalize all of the P^ by dividing each one by 
the value of S obtained in step 4 above. The P^ ' 's now 
sum to one. 

6. The values of the P-^ just obtained are used in 

rb 

the computation of P^ 114 ^ during the (n+1) th iteration, 
which proceeds exactly the same as the nth, following steps 
1 through 5 . 
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Eventually, the algorithm will converge on the solution 
of the system of balance equations.'*' The bivariate model 
balance equations were solved in the above manner. The 
average number of iterations required for convergence was 
approximately f orty-to-fif ty , with one system converging 
in seventeen iterations . The end result of the algorithm 
is a solution for { > which was used to calculate the 
measures of congestion contained in the set P (equation 
(2.3)) . 



The calculation of the elements of P is straight forward. 

Each component of the vector elements P^, P p ^ , and P B ^ is 

a standard queueing system calculation. The computation 

of each component proceeds without difficulty, once the set 

{P r , } is known, 
fb 

Define {P^} as the univariate probability distribution 
of the flight-line portion of the queueing system under 
consideration. Then 



{P f } = {P [F = f]}, f = 0,1,2, ...,M , (3.10) 



where 



The fact that the Gauss-Seidel iterative method 
described converges has been established by J.P. Lehoczky 
(private communication) . 
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P f = P[F = f] 



(3.11) 



M-f 

= E P 
b=0 



fb 



Similarly/ for the back-shop, 



P b = P [B = b], b = 0,1,2,.. . ,M , 



(3.12) 



where 



P b = P[B = b] 



M-b 

Z 

f=0 



fb 



(3.13) 



The univariate probability distribution of the number 
of aircraft in the queueing system in the aggregate 
(flight-line and back-shop) is defined as P n , where 



P n = P [N = n] , n = 0,1,2,.. . ,M , 



(3.14) 



with 



P = P[N = n] 
n 



Z P* , 
f=Q f , (n-f) 



(3.15) 



or 



P 

n 



b=0 p '»- b '- b • 



(3.16) 
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The expected number of aircraft down for flight-line 
repair at the queueing system under consideration (shop 
i) is 



M 

E [N„ . ] = £ f Pf (3.17) 

Fl f=l f 

where P^ is the value obtained from equation (3.11) . 
Similarly for the back-shop, 

M 

E[N . ] = £ b P, (3.18) 

Bl b=l b 

where P^ is the value obtained from equation (3.13). The 
expected number of aircraft down in the system as a whole 
is 

M 

E [N . ] = £ n P„ (3.19) 

i n n 

n=l 

where P n is obtained from either equation (3.15) or 
equation (3.16). 

The expected number of aircraft in the flight-line queue 
is 



M 

E [Q„ . 1 = £ (f - S.) P f . (3.20) 

Fl f=s i +i 1 f 



The expected number of aircraft in the back-shop queue is 
obtained after some head scratching. When the number of 
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aircraft that are down on the flight-line is less than or 
equal to the number of repairmen at the shop (f <_ S^) , 
f repairmen will be working at the flight-line, with the 
remaining - f men available for back-shop repairs. As 
long as the number of aircraft down for back-shop repairs, 
b, is less than or equal to s^ - f, no back-shop queue will 
exist, since there is an adequate number of repairmen 
available to handle back-shop tasks. This implies that the 
queue at the back-shop exists when b >_ - f + 1. The 

queue will contain (f + b - S^) aircraft. Thus the expected 
number in queue at the back-shop when 0 < f < is 

S i M-f 

Z Z (b+f-S.) P- 

f=0 b=S . -f+1 1 1 ' 

l 



When f _> S^, all repair personnel are tied-up with flight- 
line repairs, so that every aircraft down for back-shop 
repair is in queue, awaiting service. This yields 



M-l M-f 
Z Z 

f=s i +i b=l 



b P 



fb • 



Adding the two expressions together gives the expected 
number in the back-shop queue: 






S i M-f 

I £ (b + f-S . ) P 

f=0 b=S . +l-f 1 fb 

1 



M-l M-f 

+ Z Z b P,, (3.21) 

f=S.+l b=l 1 ‘ 
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The expected number of aircraft in queue for the system 
as a whole is similar to the back-shop queue in regard to 
derivation. When f _< the only aircraft in queue are 

those in the back-shop queue, since there are enough 
repairmen to handle flight-line repairs. Now, when f > S^, 
queues can exist at both the flight-line and the back-shop. 
The number of aircraft in queue will be (f + b - S^) . This 
yields 

S i M-f 

E[Q. ] = Z Z (b + f - S.) P-. 

1 f =0 b=S i +l-f 1 rD 

M M-f 

+ E Z (b + f - S.) P,, . (3.22) 

f=S i +l b=0 1 

The variances of the number of aircraft in the flight- 
line, back-shop, and total system are Var[Np^], Var[N B ^], 
and Var[N^], where 



Var [N pi ] 


= E[N Fi 2 ' ‘ 


- (E [N pi ] ) 2 


(3.23) 




M 2 


M 2 






= z f 2 P f 


- [ Z f P-] , 


(3.24) 




f=l 


f=l 




Var [N Bi l 


= e[n b1 2 1 - 


- (E[N B± ]) 2 


(3.25) 




M 


M 2 






= * b P b 


- [ S b P b ] , 


(3.26) 
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and 



Var [N . ] = 
1 



E[N, ] 



- (E [N ± ] ) 



(3.27) 



M 

£ 

n=l 



2 



n 



P 



n 



M 

[ l 

n=l 



n P 

n 



2 



(3.28) 



The variance of the number of aircraft in queue at the 
flight-line is 



Var[Q F .] 



E[Q Fi ] 



2 

([Q Fi ]> 



(3.29) 



M 

= £ (f-S.) 

f=S.+l 1 

l 



2 



M 2 

[ £ (f-S.) P f ] . 

f=S.+l 1 

1 (3.30) 



The variance of the number of aircraft in queue at the 
back-shop is 



Var[Q B .J = E[Q Bi 2 ] - 



S i M-f 

= £=0 b=S.+I-f ^ ^ 

1 ' 

M-l M-f 

+ £ £ b P„ 

f=S i +l b=l tD 

2 

- (E[Q Bi l) . (3.32) 



(E[Q b .]) (3.31) 

(b+f-S,) 2 P-. 
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The variance of the number of aircraft in queue at 



the system as a whole is 




2 



Var [Q i ] 



(E[Q i l) 



(3.33) 



S i M-f 
f=0 b=S.+l-f 




2 



l 



+ Z Z 

f=S i +l b=0 



M M- 



M-f 

Z (f+b-S i ) 



2 



P 



fb 



2 



(E(Q i l) . 



(3.34) 



The calculation of the expected waiting times E[W F ^], 

E [W_ . ] and E[W.] is based on the queueing relationship 

D 1 1 

"L = AW" , where L is the expected number of customers in 
the system, A is the average arrival rate, and W is the 
expected waiting time. Denoting the average arrival rates 
at the flight-line and back-shop by A Fa and A Ba , respectively, 
and the average arrival rate to the total system by A , 
the queueing relationship is restated in terms of notation 
applicable to the problem at hand. Thus, we have 





(3.35) 





(3.36) 
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and 



E [N . ] = X E [W. ] . 

_L Cl -L 



(3.37) 



Given that the state of the system is (f,b) there are 

exactly (M - f - b) aircraft that can fail into the system. 

They each have individual failure rates X . , so that the 

failure rate into the system in state (f,b) is (M - f - b) X„. . 

i? 1 

The proportion of time that the system is in state (f,b) 
is Pfb* Summing up the failure rate into the system at 
state (f,b) times the proportion of the time that the system 
is in state (f/b) over all possible states yields the average 
failure rate for the system under consideration. Thus, 



X 



Fa 



M-f M 

E E (M-f-b) X_. P (3.37) 

b=0 f=0 Fl fb 




M-f 

E 

b=0 



M 

E M P 
f =0 



fb 




M-f 

E 

b=0 



M 

E (f + b) P_ 
f=0 rD 



M-f 




[M E 
b=0 



M 

E P 
f=0 



fb 



M-f 

E 

b=0 



M 

E f P 
f=0 



fb 



M-f 

E 

b=0 



M 

E b P 
f=0 



fb 



] 



= X pi [M - E[N p .] - E[N B .]]. 



(3.38) 



X 



Fa 



= X pi [M - E [N ± ] ] 



(3.39) 
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Similarly, 



x B a = X Bi tM - EtN i 1) ' 



(3.40) 



Derivation of the overall system average arrival rate 
yields 



M-f M 

= 2 E (M - f - b) (A p . + X R .) P fb 

a b=0 f=0 Fl Bl fb 



(3.41) 



M-f M 

2 2 (M - f - b) \ . P f , 

b=0 f=0 F fb 



M-f M 

+ 2 2 (M - f - b) X . P,, 

b=0 f=0 Bl fb 



X Fa + X Ba 



(3.42) 



which is no surprise. 

The calculation of delay times is based on the queueing 
relationship Q = XD. Stated in terms of applicable notation 
the delay times are: 



E[Q Fi ] 


X Fa E *- D Fi^ ' 


(3.43) 


E ^ Q Bi^ 


X Ba E ^ D Bi^ ' 


(3.44) 
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and 



E[Q. ] = A E[D. ] . (3.45) 

X a. X 

The bivariate model is an approximate model of the 
current aircraft repair facility described in Chapter II. 

The purpose of the development of the model is to provide 
a basis of comparison for scaled-up versions such as the RAA 
and the Queen Bee structures. Two other approximate models, 
called decomposition models, were developed to assist in 
modelling the scaled-up structures. A discussion of these 
models is presented in the next chapter. 
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IV. TWO DECOMPOSITION MODELS 



He was in logic a great critic, 

Profoundly skilled in analytic. 

He could distinguish, and divide ... 

-Samuel Butler, Hudibras 

Initially conceived to provide a rough-and-ready model 
of the current repair facilities, the decomposition models 
were soon abandoned in favor of the more refined bivariate 
model described in the last chapter. However, if scaling-up 
is to occur, the current repair shops must be physically 
decomposed in some manner, as repair personnel will be taken 
from shops as currently organized and reassigned into Queen 
Bee or RAA Support Mission Bases. One scheme is to maintain 
sufficient personnel at each shop in the facility to handle 
flight-line repairs only, and reassign the remainder to an 
aggregated facility with exclusive responsibility for back- 
shop repairs. Men would be drawn from several bases for the 
aggregated back-shop. The idea is that the pooled back-shop 
resources from the bases would provide more than sufficient 
manpower to give adequate back-shop support. The excess 
manpower would then represent the manpower savings achieved 
by reorganization into the new structure. 

Part of the problem is to determine how many men to 
leave at the flight-line facilities. The bivariate model 
helps in this, since it provides a measure of flight- line 
performance at the shops represented by the vectors 
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i = 1,2,3, In addition, the bivariate model gives 
measures of back-shop performance in the vectors P* . , 
i = 1,2,3, A decomposition model fits in nicely at 
this juncture in the overall problem. With a model that 
yields measures of performance for various allocations of 
shop personnel between exclusive flight-line duties and 
exclusive back-shop duties, all that is needed is to find 
that allocation that meets present flight-line performance 
criteria, and the numbers available for the pooled back-shop 
are in hand. The back-shop can then be modeled, utilizing 
one of the simpler queueing models, to find the number of 
personnel required to provide comparable support to that 
provided by the current facility, and the manpower not needed 
represents the savings acquired by reorganization. 

Figure 6 best describes the decomposition approach. 

Again, independence between shops and exponentially distri- 
buted repair and interarrival times are assumed. A further 
assumption is that flight-line and back-shop activities within 
each shop are independent of each other. The men assigned 
to shop i are allocated so that S_. men conduct flight-line 

r 1 

repairs exclusively, with the remaining S_. . working in the 

B 1 

back-shop . 

Figure 7 is a transition diagram of the process. At 
each shop, two independent processes operate simultaneously. 
Each process is identical to that described by Gaver and 
Thompson in [Ref. 3] as the Finite Arrival Source, or 
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DECOMPOSITION MODE.L OP CUR.R.EWT FACILITY 



TRANSITION DIAGRAM OF HYPOTHETICAL 
REPAIR SHOP MODELED AS TWO INDEPENDENT 
FINITE ARRIVAL SOURCE QUEUES 
DECOMPOSITION REPAIRMAN MODEL 
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FIGURE 



Repairman problem. The state of each process can be repre- 
sented by (F ( t) ) , or (B ( t) ) as applicable. F(t) is the 
number down for repairs at the flight-line, and B(t) is the 
number down for back-shop repairs. Let 

P f (t) = P [F ( t) = f] . (4.1) 

Again, the long-run probabilities are the sought quantities, 
i.e. 



P f = lim [P (t) ] 
t-*» 

= lim P [F ( t) = f] , (4.2) 

t-*-°° 



and 



P. = lim P [B ( t) = b] . (4.3) 

t»*c© 

Balance equations for the queueing system represented 
by flight- line of Shop i are: 
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State 



Rate Process Leaves 



Rate Process Enters 



0 



1 



2 



“Sl P 0 



+ y pi ]P 1 
t (M-2) X Fi + 2v Fi ]P 2 



= MVi P o + Vi (4 ' 5 > 

= (M-U Vi P l l,+ 3 Fi P 3 (4 ' 6) 



o<f<s fl [(»-£) ^ + = (“■wi’kVi 

+ (f+1) P Fi p f+1 (4.7) 



E„. < f < M 
Fl — 



[ (M-f) Ap. + S E .y E .]P f = (M-f+1) ^p ± P f _ 1 



+ S Fi^Fi P 



f+1 



(4.8) 



M 



S Fi y Fi P 



M 






• P, 



M— 1 



(4.9) 



Replacing the term MA^Pq on the left side of the balance 
equation for state 1 with its equivalent y F ^P^ obtained from 
the balance equation for state 0, the balance equation for 
state 1 becomes 



(M-l) A pi + 



y Fi P 1 



= y Fi P 1 + 



2y Fi P 2 



(4-10) 
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Clearing terms 



(M-l) A 



Fi 



2 » Pi p 2 



(4.11) 



Replacing the expression (M-l) 1^ on the left side of the 
balance equation for state 2 with its equivalent just found 

in equation (4.11) 2 u Fi P 2 , and clearing terms, the balance 
equation for state 2 becomes 



(M-2) An • = 3)i p 

Fx 2 J ^Fi P 3 * (4-12 

This process continued down the states from state 3 to state 
M yields the following set of balance equations: 

*f P f “ y f+l P f+l ' 0 1 f 1 M_1 (4.13) 



where 



= (M-f ) A pi 


0 < f < M-l 


(4.14) 


fy Fi ' 


0 < f c S n . 
— Fl 


(4.15 


S Fi y Fi ' 


S Fi 1 f I M . 


(4.16) 


solution to this set 


of balance equations is 




f A i-1 

£ = P ° A \ ■ 


0 < f <_ M . 


(4.17) 
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Pq is obtained from the requirement that 



Thus , 



M 

Z P j 
f =0 j 



= 1 . 



(4.18) 



P„ + 



and 



M 




M 


+ Z 


p f 


= Z 


f=l 


r 


O 

II 


M 




f X . 

n 


+ Z 


p n 


l — 1 

1 II 


0 


j=i y 


i- 


M 


f X 


i + 


Z 


n - 




f=l 


L j=i 



P f = 1 ' 



- 1 , 



3 J 



j-1 

y j J 



= .1 



(4.19) 



The solution of the balance equations obtained by 
substituting expressions for X ^ and obtained from equations 
(4.14), (4.15) and (4.16) into equations (4.17) and (4.19) 
yields the follov/ing solution for the flight-line balance 
equations : 



P„ = 



S Fi 



M! 



1 + f l ± (M-f) if! 



P^ = 



Ml 



f (M-f) If 1 



P 



Fi 



i 1 


f M 

M 1 

IV 


Vi ' 


^Fi 1 


1 <M-f) !S Fi .S F1 t -S Fi 
Fi 


[Vii 

( 


\ f 






p c 


, ' 0<f i S Fi 


t 



52 



M! 



(M-f) ! S F .! S F . f S Fi 




S_. < f < M 

Fi — — 



(4.22) 



The solution for the back-shop balance equations is 
obtained employing the same method used to solve the flight- 
line balance equations. The solution is 



P 0 = 



Bi 

1 + Z 
b=l 



M! I 




b M 

+ V M! .... 


*Bi 


(M-b) lb I 


i 


b=S Bi+ l <(H>HS Bi !S Bi b - S Bi 


^Bi 



(4.23) 



M! 

P b = (M-b) !b! 



Bi 

J Bi 



° < b < s Bi 



(4.24) 



and 



P b " 



M! 



(M " b) iS Bi !S Bi b " SBi 



Bi 

J Bi 



S n . < b < M 
Bi — — 



(4.25) 



Calculation of the elements of P„. and . were made 

Fi Bi 

utilizing equations from Chapter III where appropriate. For 
the expected number of aircraft down, equations (3.17), (3.18), 

and (3.19) were used. The expected number of aircraft in 
queue at the flight-line was calculated from equation (3.20) . 
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The expected number of aircraft in the back-shop queue was 
calculated using 



E(Q 




M 

£ (b-S,) 



(4.26) 



The variance of the number of aircraft in the flight-line 
and back-shop were calculated using equations (3.24) and 
(3.26), respectively. The variance of the number of aircraft 
in flight-line queue was computed from equation (3.30). The 
variance of the back-shop queue was computed from 

M 

Var (Q_ . ) = I (b - S T ) 2 P. - (E[Q r .]) 2 . (4.27) 

Bl b=S i +l I b Bi 

Calculation of the expected waiting times and expected 
delay times were made utilizing equations (3.35), (3.36), 

(3.43) , and (3.44) with 



*Fa = *Fi [ M ' E[N Fi>] 

and 

*Ba " X Bi [ M - EIN Bi ] ] 



(4.28) 



(4.29) 



Another decomposition model assumes a Poisson arrival 
process, with constant arrival rate MX^ . or MX B ^ as applicable. 
Although not as accurate as the repairman model, it can be 
applied to systems with a large number of aircraft and 
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sufficiently manned so that the number of aircraft operational 
is kept at high levels. 

Again the problem has been reduced to the repeated 
solution of a simple model, this time the Poisson-Arrival , 
Parallel-server, Exponential Service Time System, or M/M/S 
Queue described by Gaver and Thompson in Ref. 4. The long- 
run probability distributions {p^} and (p^) can be found, 
utilizing 



P 



f 




f 



0 < f < S 



Fi 



(4.30) 



and 




f 



S 



Fi 



< f 



(4.31) 



where 




(4 32) 
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For the back-shop, use 



bO 

bi 



'mA .\ b 

Bl 
y Bi 



, 0 b _5. ^ ^ 



bO 



S i Q k-Sgi 
S Bi* S Bi 



'MA 



Bi 



Bi 



' S Bi ^ b 



(4 



where 



bO 






S Bi , 
y JL 

b=0 b! 



MA 



/ MA • \ 


b 


MA . 'l 


Bi l i \ 


Br 


4- 


Bi 


/ 1 


\ yBi 1 


T 


y Bi 1 





Bi 



Bi Bi 



MA 



1 - 



Bi 



S Bi y Bi 



-1 



// 



( 4 , 



The expected number of aircraft down for repairs at the 
flight-line and back shop are 



EtN Fi ] = 



and 



MA Fi 

— SjL + p 

%i f0 



S p . ( MA Fi 

! MX Fi\ I 1 \j S Fi y Fi 

iv) l S Fi« 

{ \ S Fi 



MA Fi \ 2 
S Fi y Fi / ; 



>, (4 



MA. 



EtN Bi ] = 



Bi 



y Bi 



MA 



Bi 



+ P. 



Bi 



bO 



^3 



MA 



Bi 



S Bi y Bi 



S Bi l 



MA 



1 - 



Bi 



>• (4 



S Bi y Bi, 



33 ) 



34 ) 



35 ) 



. 36 ) 



. 37 ) 
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The expected queue lengths are 



MX Fi 

E [Q_. ] = E [N . ] — 

Fi Fi 

Fi 



and 



EtQ Bi ] = E[N e .J 



MX Bi 

U Bi 



The expected delay times are 



EtD Fi 1 



e[q f1 ] 

HX Fi 



and 



E[D Bi ] 



E[Q Ei l 

~ a zr 



The expected waiting times are 



E[W pi ] 



E [D . ] + -=- 
Fl ^Fi 



and 



E tW Bi ] = E [D Bi ] + - 



Bi 



( 4 . 38 ) 



( 4 . 39 ) 



( 4 . 40 ) ' 



( 4 . 41 ) 



( 4 . 42 ) 



( 4 . 43 ) 
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V. COMPARISON OF MODELS 



Comparisons are odorous . 

-William Shakespeare, Much Ado About Nothing 

FORTRAN programs were written to obtain solutions for 
the bivariate and decomposition (repairman) models. Results 
for the decomposition (M/M/S) model were obtained from 
previous M/M/S queue calculations made at RAND. Data used 
in the programs was provided by RAND. It consisted of an 
arrival rate, a service rate, and the number of servers at 
each shop of a 20 shop repair facility. It was decided to 
work with only 10 shops of the 20. This decision eliminated 
unnecessary repetitive calculations which would have added 
little to the model comparisons, and still provided enough 
data so that the manpower savings resulting from a consolida- 
tion of back-shop functions could be examined. 

The data could not be directly applied to the models , 
since arrival rates and service rates were computed from all 
repair transactions at each shop, without differentiation as 
to category of repair, flight-line or back-shop. A rough 
estimate of the ratio of back-shop to flight-line repairs 
(.25) provided a basis for converting the arrival rate into 
the shop into flight-line and back-shop arrival rates. Back- 
shop and flight-line service rates were estimated to be equal, 
so no conversion was necessary. The number of aircraft used 
in the decomposition models was 25. Table II shows how the 
data was converted. 
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TABLE II 



Data Conversion Table 



Data Provided 


Symbol 


Model 






Bivariate 








and 


M/M/S 






Repairman 




Arrival Rate 


X. 

l 


X Fi = .8(^|) 


X Fi = - 8X i 






X Bi ■ - 2 4 t> 


X Bi - ‘ 2X i 


Service Rate 


y i 


y Fi = 


y Fi = y i 






y Bi = ^i 


y Bi = *i 



The division of repair personnel between flight-line and 
back-shop portions is not necessary in the bivariate model. 
But the decomposition models require that some division be 
made. For a shop with S^ personnel, there are - 1 possible 
divisions that may be made. It was decided to make the 
division so that the expected number of customers in the shop 
was minimized. To find this division, the repairman FORTRAN 
program was run for each shop, with 1 through - 1 repairmen 
working in the back-shop, and then run again with 1 through 
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- 1 repairmen working at the flight-line. The expected 
number of customers for each combination of and Sg^, such 

that S F ^ + Sg^ = , was looked at. That combination 

producing the minimum expected number of customers determined 
the division of personnel. This division of personnel was 
also used in the M/M/S decomposition model. 

The results of the calculations for the bivariate model 
and the decomposition models are presented in Tables III 
through XII and Figures 8 through 29. In some cases, 
calculations for the back-shop portions of the decomposition 
(M/M/S) model were not made. However, there is still 
adequate information presented to allow a comparison of the 
models. The bivariate model predicts the best performance, 
while the decomposition (repairman) model predicts the 
second-best performance, and the decomposition (M/M/S) model 
predicts the worst performance of the three. The performance 
predicted by decomposition (repairman) model is in most 
cases very close to the performance precicted by the bivariate 
model, with the exception of queue lengths and delay times. 

The bivariate model shows the flexibility of the queueing 
systems presently employed by the USAF, where shop personnel 
can be assigned to the flight-line or the back-shop as they 
are needed. Queueing will not occur until all personnel 
are busy in the bivariate model, whereas in a decomposition 
model, queues form when either Sp^ or Sg^ personnel are busy 
at the flight-line or back-shop of the system under considera- 
tion. The M/M/S queue bad performance is possibly caused by 
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the fact that no reduction in demand for repairmen occurs 
as the number in the system increases. 
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TABLE III 



Measures of Congestion for Shop 1 
Calculated from Three Different Models 





Calculations for 
\ x = 0.001976, p B1 = 0.448 
X pl = 0 o 00792 , ^ = 0 o 448 


Bivariate 
Model 
Si - 3 


Decomposition 
Models, Spi = 2, S B j = 1 


Measure 


1 Repairman 


M/M/S 


e(n fi ) 


0.4337 


0.4527 


0.4647 


Var(Npl ) 


0 o 4299 


0.4818 


N.C. 


E(Q F 1 ) 


0.001315 


0.01870 


0.0227 


Var(Qpi) 


0,00176 


0.02702 


N.C. 


e(w f1 ) 


2.2391 


2.3283 


2.3467 


e(d f1 ) 


0.00679 


0.0962 


0,1146 


e(n b1 ) 


0,1098 


0.1226 


0.1239 


Var(N fil ) 


0.1099 


0.1361 


N.C. 


e(Q B i ) 


0,0041 


0.0129 


0.01367 


Var(Q fil ) 


0.0020 


0.0155 


N.C. 


e(w b1 ) 


2.2711 


2.4935 


2.508 


E(D B1 ) 


0.0848 


0.2614 


0.2766 



N.C. - Not calculated 
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probability op kj aircraft down 
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TABLE IV 



Measures of Congestion for Shop 2 
Calculated from Three Different Models 





Calculations for 
A. b2 = 0.00139, u fi2 =0.272 
X f2 = 0 o 0056 , y F2 = 0.272 


Bivariate 
Model 
S 2 = 3 


Deocomposition 
Models, S F2 = 2, S b2 = 1 


Measure 


Repairman 


M/M/S 


e(n f2 ) 


0.5005 


0.5292 


0,5467 


Var(Np 2 ) 


0.4970 


0.5760 


o 

o 

o 


E(Q f2 ) 


0.00237 


0.0290 


0.03569 


Var(Q F2 ) 


0.00310 


0.0442 


o 

o 

• 


e(w f2 ) 


3.6938 


3.8894 


3.933 


e(d f2 ) 


0.0173 


0.2129 


0. 2563 


E(N b2 ) 


0.1280 


0.1488 


0 o 1 467 


Var(N B2 ) 


0.1286 


0.1636 


N.C. 


e(q b2 ) 


0.0115 


0.0176 


0.0188 


Var(Q fi2 ) 


0.0118 


0.0219 


N.C. 


e(w b2 ) 


3.773 


4.184 


4.216 


E(D b2 ) 


0.339 


0.5081 


0.5394 



N.C. - Not calculated 
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TABLE V 



Measures of Congestion for Shop 3 
Calculated from Three Different Models 





Calculations for 
Xfi 3 = 0.0123, = 0.379 

X f3 ~ 0 o 0492, P F3 = 0.379 


Bivariate 

Model 

S 3 = 11 


Decomposition 
Models, S f3 = 7, Sg 3 = 4 


Measure 


Repairman 


M/M/S 


e(n f3 ) 


2.7921 


2.8811 


3.292 


Var(N F3 ) 


2.4804 


2.5926 


o 

o 

o 


E(Qp 3 ) 


0.000006 


0.0097 


0.0469 


Var(Q F3 ) 


0.000008 


0.0176 


N.C. 


E(W f3 ) 


2.6385 


2.6475 


2,6767 


E(D f3 ) 


0.0000059 


0.0089 


0.0381 


E(N M ) 


0.6992 


0.7886 


0.8152 


Var(N B3 ) 


0,6797 


0,7691 


N.C. 


e(q b3 ) 


0.000055 


0.00154 


0.0026 


Var(Q B3 ) 


0,000067 


0.00211 


N.C. 


EC"’ B3 ) 


2,639 


2.644 


2.647 


e(d b3 ) 


0.00021 


0.0052 


0,0084 



NoC. - Not calculated 
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PROBABILITY OF N AIRCRAFT DOWN 
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TABLE VI 



Measures of Congestion for Shop 4 
Calculated from Three Different Models 





Calculations for 
* b4 = 0.00732. u b4 - 0. 
X F4 =0.02936 . li F4 = 0. 


,313 

313 


Bivariate 

Model 

s 4 = 9 


Decomposition 
Models. S F4 = 0, S B4 = 3 


Measure 


Repairman 


M/M/S 


E(N f4 ) 


2o 0991 


2.1505 


2.3683 


Var (Np 4 ) 


1 o 9320 


1.9962 


N.C. 


E (QF4^ 


0.000023 


0.0072 


0.02323 


Vat(Q F4 ) 


0.000031 


0.0123 


N.C. 


e(w f4 ) 


3.1949 


3.2056 


3.227 


E ( Dp 4 ) 


0.000035 


0.0107 


0.0316 


e(n B4 ) 


0.5235 


0.5752 


N.C. 


Var(N B4 ) 


0.5126 


0.5731 


N.C. 


E(Q B4 ) 


0.000128 


0.0040 


o 

o 

• 


Var(Q B4 ) 


0.000179 


0.0055 


• 

o 

o 

z 


e(w b4 ) 


3.1958 


3.2174 


N.C. 


e(d b4 ) 


0.00078 


0.0225 


N.C. 



N.C. - Not calculated 
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PRO&A&IL I TV OF N AIRCRAFT DOWN 
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TABLE VII 



Measures of Congestion for Shop 5 
Calculated from Three Different Models 





Calculations for 
A — 0.01 224 9 u b 5 = 
A — 0 o 0438 i ^ p 5 = 


0.373 

0.373 


Bivariate 
Model 
S 5 - 1 2 


Decomposition 
Models, Sp^ = 8 , S B 5 = 4 


1 


Measure 


Repairman 


M/M/S 


e(n f5 ) 


2.8213 


2.8943 


3.4362 


Var(N p5 ) 


2.5029 


2.5703 


N.C. 


e(q f5 ) 


0.00000084 


0.00214 


0.1655 


Var(Q F5 ) 


0.00000107 


0.0035 


N.C. 


E(W f5 ) 


2.6810 


2.6830 


2.8166 


E(D f5 ) 


0.00000079 


0.0020 


0.1356 


E(N b5 ) 


0.70245 


0.7959 


• 

O 

0 


Var(Ng^) 


0.6827 


0.7761 


N.C. 


E(Q B 5) 


0.00269 


0.0016 


N.C. 


Var(Q B 5 ) 


0.0524 


0.0022 


N.C. 


e(w b5 ) 


2.6811 


2.6864 


N.C. 


e(d b5 ) 


0.01039 


0.00541 


N.C. 



N.Co - Not calculated 
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probability OF n AIRCRAFT down 
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TABLE VIII 



Measures of Congestion for Shop 6 
Calculated from Three Different Models 





Calculations for 
A gg = 0o02256, y gg = 0o775 
A p6 = 0 o 0904 t y p5 = Oo 775 


Bivariate 

Model 

s 6 = 10 


Decomposition 
Models . Sp^ = 7 * 3 


Measure 


Repairman 


M/M/S 


E(N F6 ) 


2.545 


2c616 


3.000 


Var(N F6 ) 


2o 286 


2o3661 


N.C. 


e(q F6 ) 


0o00001843 


0 o 0051 


0.0839 


Var(Q F6 ) 


0.000025 


0.00875 


N.C. 


e(w f6 ) 


1.290 


1.2928 


1.3274 


e(d f6 ) 


0o00000935 


0.0025 


0.0371 


E(N'b6> 


0 o 635 


0.7161 


N.C. 


v ar(N B6 ) 


0.619 


0.7207 


N.C. 


e (Qb6^ 


Oo 0001 38 


0.00922 


N.C. 


v ar(Q B6 ) 


0.000182 


0.01365 


N.C. 


E(«B6> 


1.2903 


1.3072 


N.C. 


E(D B 6> 


0.00028 


0.0168 


N.C. 



N.C. - Not calculated 
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PROBABILITY OF N AIRCRAFT DOWkJ 
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PROBA&IUTY op n AIRCRAFT down 
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TABLE IX 



Measures of Congestion for Shop 7 
Calculated from Three Different Models 





Calculations for 
^gy — 0o0092 t V^y = 0 
^F 7 = 0© 0368 > ~ 0 


.298 

.298 


Bivariate 
Model 
S 7 = 6 


Decomposition 
Models, Spy = 4, Sgy = 2 


Measure 


Repairman 


M/M/S 


E(N f7 ) 


2 0 68 94 


3.1859 


4.936 


Var(Np 7 ) 


2.4881 


4.1567 


N.C. 


e(q f7 ) 


0.0257 


0.49206 


1.8490 


Var(Q F7 ) 


0.0493 


1.3335 


N.C. 


e(w f7 ) 


3.3880 


3.9687 


5.3654 


e(d f7 ) 


0.0323 


0.6130 


2.0097 


E(N b7 ) 


0.7401 


0.8476 


0.9069 


Var(N B7 ) 


0.7725 


1.0286 


N.C. 


E(Q b7 ) 


0.0756 


0.1019 


0.1351 


Var(Q B7 ) 


0.1335 


0.1869 


N.C. 


e(w b ?) 


3.7297 


3.8145 


3.943 


E(D b7 ) 


0.3810 


0.4587 


0,5872 



NoC. - Not calculated 
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PROBABILITY OF *4 AIRCRAFT DOWN 
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TABLE X 



Measures of Congestion for Shop 8 
Calculated from Three Different Models 





Calculations for 
X B8 = °- 00572 - W B8 - 

^•p8 = 0 o 0228 i ypg = 


0,472 

0,472 


Bivariate 
Model 
Sg = 7 


Decomposition 
Models. S p 8 = 5. S B g ~ 2 


Measure 


Repairman 


M/M/S 


w 

00 


1,1426 


1.1570 


1,2146 


Var(Npg) 


1.0904 


1.1087 


N o Co 


E ( Qf8 ) 


0.000013 


0.00126 


0,0027 


Var ( Q f8 ) 


0o 00001 7 


0.0018 


N.C, 


E(W F8 ) 


2.119 


2.1210 


2,1234 


E(Dpg) 


0,000025 


0.0023 


0,0048 


e(n b8 > 


0 o 2857 


0,3053 


N.C, 


Var(N B g) 


0,2825 


0.3135 


N.C. 


E(Q B 8) 


0,000051 


0.0060 


N.C. 


Var(Qgg) 


0,000069 


0.0078 


N.C, 


e(w b8 ) 


2,1192 


2,1613 


N.C. 


e(d b8 ) 


0,000378 


0,0426 


N.C. 



NoC. - Not calculated 
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PROBABILITY op n aircraft Down 
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PROBABILITY of n aircraft down 
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TABLE XI 



Measures of Congestion for Shop 9 
Calculated from Three Different Models 





Calculations for 
X - B9 = Oo 02216, y B9 = 
^ f9 = Oo 0884 , y F9 ~ 


0.333 

0,333 


Bivariate 
Model 
S 9 - 23 


Decompos ition 
Models, S p9 = 15, Sg 9 = 8 


Measure 


Repairman 


M/M/S 

A 


E(Nf9) 


4 o 93 24 


5,2444 


6,6397 


Var(N F9 ) 


3 o 9394 


4,1442 


N ,Co 


E(Qpg) 


0,000000 


0,000000 


0.0031 


Var(Qpg) 


0.000000 


0.00001 


NoCo 


e(w f9 ) 


3.003 


3,003 


3,004 


e(d f9 ) 


0,000000 


0,000000 


0.0014 


e(n b9 ) 


1.2490 


1,5599 


N o Co 


Var(N B9 ) 


1,1366 


1 . 4626 


N.C. 


e(q B9 ) 


0.000000 


0.00002 


N.C. 


Var(Q B9 ) 


0.000000 


0,00002 


N.Co 


e(w b9 ) 


3.0030 


3.003 


N.C. 


e(d b9 ) 


0,000000 


0.00003 


N.C. 



N.C«- Not calculated 
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PROBABILITY OF N AIRCRAFT DOWN 
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TABLE XII 



Measures of Congestion for Shop 10 
Calculated from Three Different Models 





Calculations for 
^B10 = 0,001 0 = 

VlO = 


0.389 

0.389 


Bivariate 

Model 

Sio - 3 


Decomposition 
Models, Spio=2» S B i \ q =\ 


Measure 


Repairman 


M/M/S 


E(N F io) 


0 o 3366 


0.3458 


0.3522 


Var(Npio) 


0o 3335 


0.3582 


N.C. 


e(Qfio) 


0o00005 


0.00865 


0.01029 


Var(Qrio) 


O 0 OOOO 6 


0.01155 


N.C. 


e(w f10 ) 


2.5745 


2.6366 


2.648 


E(D F i 0^ 


0o 0038 


0.0659 


0.0774 


e(n b io) 


0.0346 


0.0926 


0.0896 


Var(N Bl0 ) 


0.0845 


0.1003 


N.C. 


e(Qbio) 


0.000712 


0.00753 


0.00737 


v ar(Q B i 0 ) 


0.000807 


0.00374 


N.C. 


e(w b10 ) 


2.592 


2.7932 


2.801 


E(D b io) 


0.0218 


0.22759 


0.2304 



N.Co - Not calculated 
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PROBABILITY OF N AIRCRAFT DOWW 
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VI. CALCULATION OF MANPOWER SAVINGS 



When wilt thou save the people? 

Oh, God of Mercy! when? 

The people, Lord, the people! 

Not thrones and crowns, but men! 

-Ebenezer Elliot, Poetical Works 

One application of the models developed in the previous 
chapters is the determination of manpower savings attainable 
from a scaled-up facility. Specifically, the savings 
achievable if the back-shop functions of five identical bases 
were consolidated into a large back-shop facility was deter- 
mined. Each base was assumed to consist of 10 shops, each 
identical to Shops 1 through 10 used in the model comparisons 
made in Chapter V. Each base operates 25 aircraft with 
failure characteristics to those used in the model comparisons. 
The scheme was to leave crews at the 5 bases to handle only 
flight-line repairs. Back-shop repairs were made at the 
consolidated back-shop. 

The decomposition (repairman) model was used to determine 
the number of repair personnel to leave at each base for 
flight-line duties. The back-shops were manned from the pool 
of remaining maintenance personnel. The consolidated back- 
shops were modeled as M/M/S queues since the number of 
aircraft failing was large enough (125) so that the model 
predicts reasonably well. 

Each base was manned with a total of 87 repair personnel 
at the ten shops, with 58 assigned to flight-line duties and 
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29 assigned to back-shop duties (according to the repairman 
model) . For 5 bases there were 435 total men available for 
the base consolidation, with 290 assigned to flight-line 
duties, and a pool of 145 men available for assignment to 
the consolidated back-shops. 

It was decided to assign men to the consolidated back- 
shops so that the expected number of aircraft down at each 
shop was within 10% of the number predicted by the decomposi- 
tion (repairman) model, with shops manned as in Chapter V. 

For example, if the expected number of aircraft in the back- 
shop of Shop i was 1.00 before the consolidation, then 
5 bases would have, on the average, 5.00 men in the back-shops 
of Shop i. The standard for the consolidated back-shop would 
then be an expected number of aircraft down of less than 5.5. 

The arrival rate used at the consolidated back-shop of 
Shop i was 12 51^. The service rates were unchanged. The 
M/M/S queue model was run for various manning levels, until 
the expected number of aircraft in the shop dropped below 
the standard discussed above. Table XEII presents the result 
of the calculations of four measures of congestion made for 
the consolidated back-shop facility, and also shows the 
manpower savings achieved. 

A total of 89 personnel were eliminated by consolidation 
of back-shop repairs. This is roughly the former manning 
level of one base, or a savings of roughly 20% achieved by 
scaling-up the back-shop functions. Recall, that the total 
number of shops at the base from which the data was collected 
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was twenty. The number of men saved by reorganization should 
be roughly twice that predicted here, or a little less than 
200 repair personnel. 
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VII. CONCLUSIONS 



'Is there any point to which you wish to draw 
my attention?' 

'To the curious incident of the dog in the night- 
time . ' 

'The dog did nothing in the night-time.' 

'That was the curious incident, ' remarked Sherlock 
Holmes . 

-Sir Arthur Conan Doyle, The Memoirs of Sherlock Holmes 

The decomposition (repairman) model was in close 
agreement with the bivariate model, except with respect 
to queue lengths and delay times. It is quite suitable for 
use in problems such as base consolidation studies, where 
savings of repair personnel are sought. The decomposition 
(M/M/S) model provides conservative predictions, which makes 
it more suitable for studies made for a risk-averse decision 
maker. 

If performance predictions are desired for repair 
configurations similar to that described in Chapter II, 
where the effects of sustained changes in tempo of operations 
or the effect on system performance of manpower reduction or 
augmentation is desired, then the bivariate model is more 
appropriate than either of the two decomposition models. 
Furthermore, although all three models are over-simplified, 
the balance equations describing more complex situations can 
easily be written down. Then they may also be solved using 
Gauss -Seidel iteration. Thus techniques to replace complex 
simulation, at least in some instances, are at hand. Such 
work is under continuing development. 
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